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Cosmic microwave background power spectrum estimation with the 
destriping technique 
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ABSTRACT 

Extraction of the CMB (Cosmic Microwave Background) angular power spectrum is a chal- 
lenging task for current and future CMB experiments due to the large data sets involved. Here 
we describe an implementati on of M ASTER (Monte carlo Apodised Spherical Transform 
EstimatoR) described in lHivon et al] 12002) which exploits the destriping technique as a map- 
making method. In this method a noise estimate based on destriped noise-only MC (Monte 
Carlo) simulations is subtracted from the pseudo angular power spectrum. As a working case 
we use realistic simulations of the PLANCK LFI (Low Frequency Instrument). We found that 
the effect of destriping on a pure sky signal is minimal and requires no correction. Instead 
we found an effect related to the distribution of detector pointings, which affects the high-£ 
part of the power spectrum. We correct for this by subtracting a "signal bias" estimated by 
MC simulations. We also give analytical estimates for this signal bias. Our method is fast and 
accurate enough (the estimator is un-biased and errors are close to theoretical expectations for 
maximal accuracy) to estimate the CMB angular power spectra for current and future CMB 
space missions. This study is related to PLANCK LFI activities. 

Key words: methods: data analysis - cosmology: cosmic microwave background. 
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1 INTRODUCTION 

In the favoured model of structure formation driven by inflation, 
primordial fluctuations are expected to be Gaussian distributed. In 
this case all the statistical information encoded into CMB (Cosmic 
Microwave Background) anisotropies is completely described by 
their angular power spectrum Ci. The main issue in this paper is the 
extraction of the Ci values starting from an observed CMB map. In 
recent years a maximum likelihood approach to the problem has 
been developed (Gorski et al. 1996) and successfully applied to the 
COBE - DMR data l Benn ett et al]|l99rj) . This process, however, 
involves a number of operations scaling as ~ N pix , where N p i x is 
the number of pixels that cover the sky. 

Recent ground-based and balloon-borne experiments have im- 
proved our knowledge of the CMB anisotropy due to the better an- 
gular resolution and higher sensitivity of these experiments. Due to 
the sizes of their data sets they have already posed challenges in 
extracting an unbiased estimate of Ci. This will become more de- 
mandi ng with curren t and future sat ellite ex periments like WMAP 
iBennett et al"ll2003l) and PLANCK frauberlbood) which produce 
maps with iV p i x ~ few x 10 6 . Therefore, brute force maximum 
likelihood estimation of the CMB angular power spectrum is not 
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feasible without approximations or assumptions specific to partic- 
ular experimental or observing strategies lOh. Spergel & Hinshawl 
ll99StlWandelt & Hansenl2003l) . 

New methods are under stud y and some of the m have al- 
ready produced interesting results. iHivon et al] 12002) introduced 
the MASTER technique used to construct the B OOMERanG an- 
gular power spectrum in Netterfield et al. 1 2002). An extension of 
MASTER is being prepared to treat PLA/VCAMike data sets. The 
core of the MASTER technique is i) a high-pass filter applied to 
the TOD (Time Ordered Data) as part of a naive map-making pro- 
cess in order to reduce instrumental non-white noise, and ii) MC 
(Monte Carlo) simulations of pure signal, pure noise, and signal 
plus noise to correct for data filtering and instrumental noise, and 
for estimating the error bars of Ct. The approach has been shown 
to be unbiased and it gives error bars close to theoretical ones when 
cosmi c variance a nd instrumental noise properties are concerned 
faivon et all2002l) . 

iBalbi et alJ fc002l) follow a MASTER-like approach but in- 
stead of using a high-pass filter they consider an IGLS (Iterative 
Generalized Least Square) map-making algorithm already devel- 
oped to treat PLANCK data jNatoli et all200ll) . Furthermore, they 
estimate noise properties needed for proper MC simulations di- 
rectly from the data INatoli et al. 2002). The extracted Ci are unbi- 
ased and have close to optimal error bars. 
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In this p a per w e present an approach similar to that 
of iHivon et alJ 120021) but exploiting the destriping algorithm 
I Buris ana et alll997UDelabromlldl99&lMainoll999T:lMaino et alJ 
1 19991) in the map-making process. Although destriping algorithms 
are considered approximations of proper IGLS map-making in the 
sense that they will not necessarily produce the minimum variance 
map, they are able to provide estimates of various systematic ef- 
fects, remove drifts f rom TOD and return cleaned TOD (e.g., see 
iMennella et al J200^ for an application with periodic fluctuations). 
Furthermore, destriping ma kes no assum ptions on the beam shape. 
It has been demonstrated iMaino et alll2002h that the instrument 
noise is the driver of the destriping performance regardless of the 
beam shape (including sidelobes). 

A generalized maximum likelihood approach to the destriping 
method has been implemented l Keiha neii et all20 03) which is able 
to fit different sets of base functions (in addition to the simple con- 
stant baseline) and, in principle, could better remove the contribu- 
tions of different systematic effects from the TODs. All these prop- 
erties make destriping attractive. Additionally, it is fast and needs 
no prior information on the instrumental noise. Note, however, that 
as we combine destriping with the MASTER approach for the Ct 
estimation, this part utilizes information on noise properties. Both 
iterative and non-iterative methods to estimate the noise character- 
istics directly from the data have been proposed iDore et all200lt 
iNatoli et all2002l). These methods have been applied in the Ci es- 
timation (e.g. lNetterfield et all200llBalbi et all2 002). 

This paper is organised as follows. In Section|2|we describe 
the destriping technique. In Section[3]we review the MASTER ap- 
proach for the extraction of Ci from a map. We have then applied 
our combination of destriping with the MASTER approach to sim- 
ulations of one PLANCK LFI 100 GHz detector. The applied scan- 
ning strategy and sky coverages are described in Section|4| In Sec- 
tion [5] we explain our findings concerning the filter function and 
signal bias that are possible means to model the effects caused by 
our map-making method on the spectrum estimates. The analysis 
pipeline and the CPU times required to run it are explained in Sec- 
tion[6] SectionQpresents the simulations and the simulation results. 
We draw our conclusions in Section|8| In AppendixlAlwe describe 
in more details the signal bias approach for modelling the effect of 
the distribution of detector pointings in output map pixels. 



2 DESTRIPING TECHNIQUE 

The destriping technique for map mak ing has be en derived from 
the COBRAS/SAMBA Phase-A s tudy iBersanelli et alll996h and 
was originally implemented by Burigana et al . (1997) and by 
iMainoet alJ J 19991) . (See also lDelabrouillelll998l ) Developed for 
the PLANCK satellite it makes use of the fact that PLANCK is a 
spinning spacecraft. Detector beams are drawing almost great cir- 
cles on the sky. Each scan circle is observed 60 times before the 
spin axis is repointed. In order to reduce the level of instrumental 
noise, the signal can be averaged over these 60 scan circles. Al- 
though the destriping technique has been developed in this frame- 
work, it could in principle be applied to any experiment with an 
observing strategy with repeated measurements of the same regions 

of the sky. 

As pointed out bv ljanssen et aljil996l) . the effect of the instru- 
mental noise, in particular 1// noise, on the average scan circle sig- 
nal can be approximated by a uniform offset or "baseline". The key 
problem in destriping is to find the magnitudes of the baselines. The 
destriping technique uses the redundancy of the observing strategy, 



i.e. it considers all intersections (crossing points) between the scan 
circles, to obtain these magnitudes. We identify a crossing point 
when two points from different scan circles fall inside the same 
pixel in the sky. The small scale structure in the signal will cause 
destriping to insert artefacts at some level when the pointings of the 
observing beam do not fall exactly on each other within the pixels 
of the crossing point search. We will come back to this point later 
in Section |5] In this s tudy we utilize the HEALPix 1 pixelisation 
scheme (Gorski, Hivon & Wandelt 1999). Its pixel dimension is set 
by the A^ide resolution parameter. A map of the full sky contains 
12iV s 2 idc pixels. 

After the baseline magnitudes have been recovered they can be 
used to subtract an estimate of the instrumental 1// noise from the 
observed TOD. The final sky map is built from this cleaned TOD 
by simply coadding (averaging) the observations falling in each sky 
pixel. Some residual 1// noise will remain after the subtraction of 
the baselines. However, for knee frequencies ft ^ 0.4 Hz the uni- 
form baseline approximation works well and the level of residual 
noise stripes in th e output map is well below the white noise level 
iMaino et alJl2002l) . The knee frequency considered in this study 
(,fk =0.1 Hz) lies well below the above limit. 



3 THE MASTER APPROACH 

The CMB temperature field observed over the whole sky can be 
expanded in spherical harmonics Yi m , 

AT(n) = ^a (m y <m (n), (1) 

ira 

where the Yg m form a complete orthogonal set of functions over the 
sky. Here n is a unit vector pointing to a location on the sky. If CMB 
temperature fluctuations are assumed to be Gaussian distributed, 
the expansion coefficients ai m are Gaussian distributed complex 
random variables with expectation value 

{dim) = , (2) 

and correlation 

(a£ m a^/ m /) = S u /5 mm iCi h , (3) 

where C| h is the underlying "theoretical" angular power spectrum 
of the CMB. 

In the case of full sky observations an unbiased estimator of 
Cf 1 is given by 

i 

m = -t 

The coefficients Ci are distributed with mean equal to the theo- 
retical value Cf 1 and standard deviation, "cosmic variance", given 
by ACe = y / 2/(2£+ l)Cf (e.g. lKnoxlll995l) . Here v = 21 + 1 
is the number of degrees of freedom corresponding to a given mul- 
tipole number £. 

From a pixelised map of the full sky the coefficients aim are 
estimated as a sum over pixels 

aim = rip ^2 w p T p Y e*m (n P ) , (5) 

V 

1 http://www.eso.org/science/healpix 
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where Sl p is the pixel area (same for all pixels in HEALPix), T p 
the pixel temperature, n p the unit vector pointing to the centre of 
the pixel, and w p the pixel weight. Throughout this study the pixel 
weights equal one for all observed pixels. 

Including the effect of instrumental noise and a symmetric an- 
tenna beam, the expectation value of the "observed" angular power 
spectrum Ce is (by inserting ae m to Eq. 0) 



{Ce) = Blcf + {N t ) , 



and variance 



2 - [BjcT + m 2 



2£ + 



(6) 



(7) 



where Be is the shape of the beam in I space, and {Nt ) is the expec- 
tation value of the angular power spectrum of the noise projected 
onto the sky. 

Thus an unbiased estimate of Cf 1 is given by 



Ce = B e 



(C t - (N t ) 



which has the expectation value {Ce} = Cf 1 and variance 



2£ + l 



Bf 



(8) 



(9) 



The situation in a real experiment is more complicated since 
the entire sky may not be available for the observation. This may 
happen due to partial sky coverage of the scanning, removal of cor- 
rupted data or due to cutting away the galactic plane region. The 
pixels are therefore weighted by a window function W(n), that 
states which pixels are observed and which are not, in a top-hat 
fashion (used here) or using a more complex apodization scheme. 
On the cut sky the Ye m are no more an orthogonal set of functions 
and a correlation between the I modes is introduced. However, it 
is still possible to determine the expansion coefficients ae m over 
the cut sky, by summing over the observed pixels in Eq. The 
pseudo power spectrum Ce can be obtained by inserting the co- 
efficients to Eq. J4j. The expectation value of the pseudo power 
spectru m can be r elated to the true angular power spectrum C\ of 
the sky jffivon et all2002l) 



{Ce)=J2 Mu ' C - 



(10) 



where the kernel matrix M e e> describes the mode-mode coupling 
due to the cut sky. For the top-hat window function this matrix de- 
pends only on the geometry of the cut applied to the data. 

In our method the temperature of a pixel of the output map 
is an average of all cleaned (baselines subtracted) observations 
falling in that pixel. The pointings of the observations are dis- 
tributed across the pixel area. The distribution of the observations 
is likely to cause some distortion in the observed CMB angular 
power spectrum, that may need to be corrected in the power spec- 
trum estimation stage. In this study the signal TOD is generated by 
scanning a pixelized sky map realization. The spread of the point- 
ings is modelled by picking the temperatures from a map (the sky 
map) with a higher resolution than the resolution of the final output 
map. For this study the pointings of the scan circles between suc- 
cessive satellite spin axis repointings are assumed to fall exactly 
on top of each other. In real experiment this will not happen due 
to various non-idealities in the satellite motion (e.g. nutation of the 
spin axis pointing, Ivan Leeuwen et af]l2002l) . The way we model 
the pointing distribution may exaggerate the distortion of the esti- 
mated CMB angular power spectrum. 



In the original MASTER implementation Jfflvon et alj|2002l) 
a high-pass filter was applied to the TOD to reduce low frequency 
noise. The effect of this filtering on the pseudo angular power spec- 
trum was modelled by a filter transfer function Fe, so that we would 
have an approximate relation 



{Ce} = V M u ,F v Bi,Cf + {N e } . 



(11) 



The filter function Fe is to be determined by signal-only MC simu- 
lations. 

In our approach the low-frequency noise is removed by de- 
striping, so there is no need for a high-pass filter. However, our 
map-making method may still have some effect on the signal, that 
needs to be corrected for. In Section [5] we study whether this ef- 
fect could also be modelled with a filter function in the way of Eq. 
<1 It . Then the Ce is the pseudo power spectrum obtained from the 
destriped map, B| is a function including the beam response and 
the effect of finite output map pixel size, and {Nt) is the average 
of the noise angular power spectrum (on the cut sky, destriped). In 
Section|5|we introduce an alternative ("signal bias") which appears 
more suitable for our case. 

In a real CMB experiment a single noisy realization of the 
pseudo power spectrum Ce is available for the estimation of the 
true angular power spectrum C\ h of the sky. In this case Eq. <1 It 
can be inverted to obtain a formula for the estimate of Cj , 



C, 



(12) 



in the filter function approach. Using the cut sky introduces corre- 
lations between the Ce- 



3.1 Coupling Kernel 

The mode-mode coupling kernel M lt i is completely determined 
once the window function VK(n) is specified. iHivon et alJ ! 2002) 
have demonstrated that the kernel depends only on the angular 
power spectrum We of the sky window function applied to the data, 



4tt 



ta 




(13) 



where 



is the Wigner 3— j symbol. There- 



at ta h 
mi 7712 m.3 

fore, once the window function W^(n) applied to the data is given 
it is a straightforward task to compute the corresponding coupling 



simulations involved in our technique and needs to be computed 
only once. 



3.2 Modelling of the Instrument Noise 

A vital part of the MASTER approach is the estimation of the an- 
gular power spectrum of the instrument noise {Ne) (see Eqs. <l 1 1 
and <12» . If the estimate of {Ne} is incorrect a biased Ce estimate 
will be obtained. 

Several methods have been developed to estimate noise prop- 
erties directly fr om the data by making us e of both it erative 
iDore et aj200ll) and non-iterative algorithms iNatoli et all2002h . 
If noise properties can be derived from the data and reasonably 
modelled, it is possible to estimate (Ne) by making use of MC 
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Figure 1. Angular power spectrum of the instrument noise. Top panel: 
Solid horizontal line is the theoretical white noise spectrum. The upper 
curve is the noise spectrum including 1/f noise for a single noise real- 
ization and no destriping (naive coaddition from TOD to map). The middle 
curve is the power spectrum after destriping. Middle panel: The section 
indicated by the box (top panel) is shown here. Six blobs are evident in the 
no-destriping curve. Bottom panel: The top curve is the mean power spec- 
trum (of 100 MC realizations) after destriping. The solid horizontal line is 
the theoretical white noise spectrum. All noise spectra refer to the mod- 
elled full set of 24 PLANCK LFI 100 GHz radiometers and nearly full sky 
coverage. 



realizations. In this approach a set of Nj£q pure-noise TOD real- 
izations are generated. They are then processed with destriping and 
coadded into maps with galactic cut applied (if relevant). There- 
after, the angular power spectra N( will be extracted. Finally we 
can construct the MC average that can be used as an estimate of 
(Ni) in Eq. fT2l . 

In this study we do not address the question of estimating the 
noise from the data. Thus, whereas in a more realistic simulation 
the noise power spectrum to be used for the noise MC simula- 
tions would be estimated from the data as discussed above, we 
have used the same underlying noise power spectrum P(f) for 
the noise MC simulations as was used for the signal+noise TOD 
from which we are estimating the Ct. This was a mixture of white 
and 1/f (slope = — 1) noise with knee frequency fu = 0.1 Hz. 



The nominal rms (root-mean-square) level per integration time 
(t = l//sam P ic = 1/108.3 s) is ~ 4.8 mK (thermodynamic tem- 
perature scale) for one 100 GHz LFI radiometer. The full set of 
24 radiometers at 100 GHz was modelled by reducing this num- 
ber by a factor of v24. We used an SDE (Stochastic Differential 
Equation) 2 algorithm to generate the TODs of instrumental noise. 
It is unrealistic to extend the 1/f noise spectrum behaviour to 
very low frequencies. Therefore we have set a minimum frequency 
/mm = 4 x 10~ 6 Hz (~ (70 h)- 1 ) below which the noise spec- 
trum becomes flat. We used AT^c = 100 noise TOD realizations 
to produce the estimate of (Nt). According to lBalbi et alJ i20o"3) . 
already 22 MC realizations produce a good estimate. 

If the instrument noise appears to be non-stationary, SDE can 
be applied to generate the overall TOD stream in a piecewise man- 
ner where each piece is assumed to be stationary with its own spe- 
cific set of noise model parameter values. The accuracy of the de- 
termination of the baseline magnitudes in the destriping is not de- 
pendent on whether the noise is stationary throughout the TOD or 
only piecewise stationary. 

Fig.Qshows the angular power spectrum of the noise map (no 
destriping) compared with the angular power spectrum after the de- 
striping. Two things are evident. Firstly, without destriping there is 
a distinct blob structure which is caused by the geometry of the 
scanning. The amplitudes of the blobs depend on the value of fk- 
The second thing is the significant capacity of the destriping in re- 
moving the blobs and making the spectrum approach the spectrum 
of the white noise. The spectrum deviates from the spectrum of the 
white noise at low I, This is an indication that some residual 1/f 
noise remains after the baseline subtraction. 



4 SCANNING STRATEGY AND WINDOW FUNCTION 

The order in which the samples of the TOD hit the pixels of the map 
is determined by the scanning strategy. Therefore any artefacts pro- 
duced in the maps due to the map-making method can be expected 
to be affected by the scanning strategy. In this study we used the 
nominal PLANCK scanning strategy. 

In this strategy the spin axis of the satellite follows the ecliptic 
plane. The spin axis is kept anti-solar by repointing it by 2.5 ar- 
cmin every hour. The spacecraft rotates around the spin axis at a 
nominal rate of 1 rpm. During one hour PLANCK scans the same 
circle on the sky 60 times. The scanning pattern used in this study 
corresponds to the 100 GHz LFI detector number 10 (location 
(6,4)) = (3? 737, 126? 228) in the focal plane image). The angle 
between the satellite spin axis and the optical axis of the telescope 
is 85°. We assumed no spin axis precession. The antenna beam pat- 
tern has been modelled by a symmetric Gaussian beam with nomi- 
nal resolution of FWHM = 10 arcmin (Full Width Half Maximum). 

Since we assume idealized satellite motion, where the 60 cir- 
cles between repointings fall exactly on each other, sample by sam- 
ple, we can immediately average these circles into a single ring. 
Thus our simulated TOD consists of 5040 scanning rings, corre- 
sponding to 7 months of measurement time, with 6498 samples on 
each ring, corresponding to a sampling frequency of /sample = 
108.3 Hz. 

The sky coverage of the applied scanning strategy was / s k y = 



2 SDE is one of the two methods in the PLANCK Level S pipeline for 
producing simulated instrument noise. The method was implemented by 
B. Wandelt and K. Gorski and modified by E. Keihanen. 
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Figure 2. The maps with numbers of hits per pixel for the nearly full sky 
(top) and for the galactic cut (bottom) cases. The maps are in ecliptic coor- 
dinate system. Note that the scale is log 1 ( ) (ni 1 it), where n^it is the number 
of hits in a pixel. 

0.985, leaving areas close to the poles unobserved. This case is 
called "nearly full sky" in this paper. In addition, we considered a 
case called "galactic cut", where the galactic region (|6| ^ 20°) 
is cut out from the maps, leaving f s ^ y = 0.646. The numbers of 
hits in the output map (N s ide = 512) are shown in Fig.|2| The polar 
regions with no hits and the galactic cut (lower map) are clearly vis- 
ible. The average number of hits per pixel was 10.6 in the observed 
part of the sky ("nearly full sky" case). 

The window functions W^(n) applied in the nearly full sky and 
in the galactic cut cases were top hat window maps, where each un- 
observed or cut out pixel was assigned a value zero and for the rest 
of the pixels value one was assigned. The resolution of the window 
map was the same as the resolution of our output map. The angular 
power spectra (We) of these window maps were determined with 
the ANAFAST code of the HEALPix package. With the aid of the 
angular spectra the kernel matrices were produced (see Eq. <13» . 



5 FILTER FUNCTION AND SIGNAL BIAS 

Often some kind of filtering is applied to the data to reduce the 
effect of low-frequency noise. The filtering is applied to the TOD 
or to the maps. In our case the TOD is "filtered" in a sense that 
a baseline for each scan circle is subtracted. The magnitudes of 
the baselines are mainly determined by the 1 // noise. However, as 
described in Section|2|the baseline magnitudes can be influenced by 
the CMB signal as well and therefore the CMB signal is affected by 
the baseline subtraction. This effect on the angular power spectrum 
has to be accounted for. The filter function Fe was introduced in Eq. 
<11> as one possible means to model this effect. The filter function 
will also include the power spectrum distortion caused by any other 
effects inherent in the map making method. 

We determined the filter function Fe in the following way. 
Starting from a theoretical input CMB angular power spectrum C| 
we produced = 450 pure signal (i.e. no noise) full sky map 

realizations with A^idc = 1024 by using the SYNFAST code of 
the HEALPix package. These maps were then observed using the 
nominal PLANCK scanning strategy 3 (see Section|4}- 

The resulting TODs were destriped and coadded into maps. 
The angular resolution of the crossing point search and the reso- 
lution of the output maps were both 7 arcmin (corresponding to 
A'sidc = 512). The pseudo angular power spectra (Ce) were then 
extracted from the output maps. The Ce were also determined from 
maps with no baselines removed (no destriping). By comparing the 
filter functions with and without destriping we can distinguish from 
other effects those filter function contributions which are due to the 
baseline removal. 

The average of these MC realizations, (CV)mc, was then used 
in place of (Ce) in Eq. 1 1 1 i and Fe was solved from this equation 
(with (Ne) = 0). 

In order to find out whether the filter function Fe depends 

3 Simulation Software is part of the Level S of the PLANCK 
DPCs and is available for the PLANCK collaboration at 
http://planck.mpa-garching.mpg.de 
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Figure 3. The input CMB power spectra Cf 1 of the ACDM model 
(solid curve) and the OCDM model (dotted curve). The ensemble average 

fs) 

(Ny^Q = 450) of the deconvolved pseudo power spectrum from signal- 
only simulations of the ACDM model is shown as well (dashed curve). The 
spectrum was derived from the destriped maps and the sky coverage was 
nearly full sky. The filter function is the ratio of the deconvolved power 
spectrum (dashed curve) and the input power spectrum (solid curve). 



0.95 



DESTRIPING, NEARLY FULL. SKY 


/ / 


ACDM (BLACK) 




OCDM (GREY) 


/ / 
/ / 




/ / 




/ / 





200 400 600 800 1000 1200 1400 



Figure 4. Filter functions Fe for two widely different input CMB power 
spectra C| . The sky coverage is nearly full sky. The resulting filter func- 
tions after the baselines have been removed from the TODs (destriping) are 

(s) 

shown. The filter functions are results of N^ c = 450 MC realizations. 
The filter functions of destriping and no destriping are nearly identical. 

on the sky signal we considered two widely different cosmolog- 
ical models: a ACDM (cosmological constant + Cold Dark Mat- 
ter) model and an open CDM model (OCDM). Their power spec- 
tra C\ h are depicted in Fig. [3] They were computed from the re- 
specti ve cosmological mod els by using the CMBFAST code (see 
Selia k & Zaldarriagalfl99r3 and references therein). As an exam- 
ple, the ensemble average of a deconvolved pseudo power spectrum 
J2e> MwiCe^Aic/Be in the case of nearly full sky and ACDM 
is shown in Fig.|5| The filter function is a ratio of this spectrum and 
the theoretical power spectrum C\ h (compare to Eq. <12> ). 

The resulting filter functions with destriping and for the nearly 
full sky are shown in Fig. |4| (The galactic cut case looks similar, 
just somewhat noisier.) In order to see in more detail the effect of 
destriping the ratio between the filter functions with and without 
destriping were derived. They are shown in Fig. [5] A number of 
things are evident (see Figs.|4]and[5j: 

(1) The filter function stays close to one up to I ~ 700. The de- 
4 http://www.cmbfast.org 
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viations from one are completely dominated by remaining cosmic 
variance (reduced here by (-^mc) -1 ^ 2 ^ ue to avera g m g over ^mc 
realizations). Thus we conclude that the effect of our map-making 
method on the angular power spectrum is very small for I < 700, 
and we could just use 1 for the filter function, i.e., no filter function 
is needed for this part of the spectrum. 

(2) For / > 700 the filter function deviates from 1 significantly. 
This high-£ tail of the filter function is strongly dependent on C\ h . 

(3) The filter functions with and without destriping are nearly 
identical. This shows that the high-^ tail is caused by some other 
effects than the baseline removal. The largest relative impact due 
to baseline removal is around 0.2 per cent (see Fig. [3}- The ratio 
contains blob structure due to low-level stripes introduced by the 
baseline removal (cf. Fig.0. 

Before drawing further conclusions we need to investigate 
which part of the high-£ excess power of the deconvolved pseudo 
spectrum (see Fig. [3) is caused by numerical errors introduced by 
the procedures we use for generating sky maps, deriving pseudo 
spectra from the maps and deconvolving the pseudo spectra with 
beam, pixel weight and mode-mode coupling kernel matrix (M u i). 
For that purpose we produced 100 realizations of the full sky CMB 
maps from the ACDM input power spectrum, applied our "galac- 
tic cut" window function, and determined the pseudo spectra of the 
resulting maps. The mean of these pseudo spectra was then decon- 
volved with the beam, with the pixel weight (provided by HEALPix 
package) and with the kernel matrix. The resulting deconvolved 
spectrum had a close resemblance to the theoretical input spectrum 
Cj h . Their ratio is depicted in Fig.|6|(black curve). The ratio stays 
close to one at all multipoles considered, showing that the proce- 
dures we use for the generation of sky maps and pseudo spectra and 
for the deconvolution are not responsible for the high-^ tail of the 
filter function (shown in Fig.|4j- 

The dependence of the filter function on Cf 1 (see Fig. [4j 
shows that it is not a good way to model the effect of our map- 
making method. That the effect appears at high t suggests that it 
might be related to the pixel size used. Therefore we determined the 
filter functions also for map resolutions ~N<^& C = 2048 (sky map) 
and iV s idc = 1024 (output map). Results are shown in Fig.Q a 
high-i? tail still exists, and is in fact somewhat higher than in the ear- 
lier case of iVside = 512 (output map). Thus the filter function tail 
cannot be removed simply by increasing the map resolution. (Note 
that the N B id e = 512 pixel size is about 7 arcmin and should thus 
resolve the power spectrum up to I ~ 180° /7 arcmin ~ 1500.) 

The binning o f the detecto r pointi ngs in pixels has been treated 
in time domain by iDore et alJ 1200 ll) . In a noise-free experiment 
they considered a residual between an observed TOD and another 
TOD produced by scanning a pixelized map that is naively coad- 
ded from the observed TOD. This residual was called pixelization 
noise. 

We have determined (Appendix^) that the main effect caus- 
ing the high-i? tail of the filter function is related to the distribu- 
tion of the detector pointings on the sky and its relation to the pix- 
elization. Comparing the filter functions of the ACDM and OCDM 
cases to the corresponding C\ one notices that they appear to be 
inversely related. Since the filter function is essentially a ratio be- 
tween (Ci) and C\ , this suggests that the dominant contribution 
would be better modelled as an offset St, which we call signal bias, 
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Figure 5. The ratio between filter functions with and without destriping. 
Nearly full sky coverage and ACDM cosmological model were applied. 
Note that different panels show the ratio at different ranges of I. 
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Figure 6. CMB power spectra were produced by first generating 100 real- 
izations of the full sky, then cutting out part of the observations (galactic 
cut case, see Fig. 0, then taking the mean of the pseudo power spectra 
derived from the partial sky maps and finally deconvolving the mean with 
beam, pixel weight and kernel matrix. The ratio between the deconvolved 
spectrum and the theoretical input spectrum Ci" is shown (black curve). 
The ratio without deconvolution with the kernel matrix is demonstrated 
as well (grey curve). The cosmological model was ACDM. All maps had 
N Bide = 512. 



function. An equation similar to Eq. J 121 can be written for the 



(C e )=y M w Blct + St, + {N e ). 



(14) 



estimate of C\ , 

^ ( 
Ct = ^ 



(15) 



It can be estimated by signal-only MC simulations just like the filter 



The average pseudo power spectrum of the MC realizations, 



CMB Power Spectrum Estimation with Destriping 7 




200 400 600 800 1000 1200 1400 

/ 



Figure 7. Same as Fig. |4] but now the sky map realizations (scanned for 
TOD) and the output map had higher resolutions (N s ^ c = 2048 and 

fs) 

Aside = 1024, respectively). The filter function is a r esultof N$ c = 450 
MC realizations. 

(Ce)MC, can be used in place of (Ce } in Eq. J 141 and S^ 10 can be 
solved from this equation, with (Ni) = 0. 

Further details of the signal bias are derived in Appendix^ It 
is shown there that the signal bias estimates are rather similar even 
though the cosmological models are widely different. 

The signal bias is reduced when data from several detectors 
are combined, since that increases the density of detector point- 
ings. We find in Appendix |A] that the signal bias in the region of 
high £ is nearly inversely proportional to the number of detectors. 
In reality, the detector pointings from successive scan circles also 
do not fall exactly on top of each other. Therefore, in the case of 
the real PLANCK experiment this signal bias may be small enough 
(when compared to other systematic effects) to be ignored. 

In our idealized simulation, which utilizes only pointings of a 
single detector and assumes exact overlap of successive scan cir- 
cles within the same pointing period, the signal bias is larger, and 
shows up prominently in our results. For comparing the signal bias 
to other map-making artefacts we carried out destriping on a noise- 
less TOD containing CMB and galactic foreground signals. In this 
case the signal bias was significantly larger (close to two orders of 
magnitude larger at high £) than the residual power spectrum due 
to the subtracted baselines. 

We show in Appendix [S] that the signal bias can be reliably 
estimated by MC simulation and thus removed. We have done this 
in what follows. 



6 ANALYSIS PIPELINE AND COMPUTING TIME 

The starting point for our analysis is a TOD with CMB+noise (the 
"original TOD"), and a noise power spectrum P(f). Our analysis 
pipeline is composed of the following steps: 

(i) A total of N^c pure noise TODs are produced, using P(f)- 
They are destriped and maps are coadded from the TODs from 
which the baselines have been removed. The angular power spec- 
tra of the output maps are produced with the ANAFAST code. By 
averaging the obtained spectra an estimate for (Ne) is obtained. 

(ii) From the original TOD a destriped map is produced, from 
which a power spectrum estimate C' t is obtained, using Eq. |15} 
with Si = 0. 

(s) 

(iii) Random generation of N^q pure CMB sky realizations 
with the SYNFAST code. The estimate C' t is used as an input power 



spectrum. In order to produce a TOD, the sky is observed according 
to the selected scanning strategy. TODs are destriped and projected 
back onto the sky. The maps are input to ANAFAST for the deter- 
mination of their angular power spectra. These spectra are used for 
determining the signal bias S^ 10 . 

(iv) The power spectrum estimate Ce is obtained, using Eq. G3 
withSf . 

At this point we have our final power spectrum estimate Ce- 
To obtain an estimate for the error bars and the covariance matrix, 
further simulations are required. In the case of a real experiment, 
these simulations would be based on the Ci just obtained. However, 
in this study we have also available the Cf 1 used to produce the 
original TOD. The emphasis of this study is to reveal the actual 
error bars of our method, rather than to simulate the estimation of 
these error bars in a real experiment. Thus we use this C} h as the 
input spectrum to generate an ensemble of MC estimates Ce . 

(v) A total of A^q 11 ' randomly generated "experimental" TODs 
(i.e. with sky signal and noise) are produced. They are destriped 
and coadded into a map. Thereafter the pseudo power spectra Ce 
are determined. 

(vi) The estimated angular power spectra Ce are computed using 
Eq. d 1 5 i and binned. 

(vii) The ±lcr error bars on each Ce and bin are evaluated as std 
(standard deviation) calculated over the MC realizations. 

In our simulations we had N^l = 100, N$ c = 450, and 
Af^J 11 ' = 450. The relatively large number of signal+noise MC 
cycles was required to study the statistics of the C'e estimates. The 
relative rms accuracy of the ±1<t error bars is (2A^J n ') _1 / 2 . 
With Aj^J 11 ' = 450 this implies 3.3 per cent rms accuracy. Ad- 
ditionally, this number of MC realizations is high enough to reveal 
some features of the covariance matrix C«/ = {(Ce — (CV))(CV — 
(Cey ) )) (Section^. However, to obtain estimates for individual off- 
diagonal elements of Cee' > °r even of the binned covariance matrix 
CbV = {(Cb — (Cb))(Cb> — (Ci,/))) would have required an even 
larger N^, n ' > . For A^J 11 ' = 450 the off-diagonal elements were 
still dominated by MC noise for bins of A£ = 10. 

As for the original MASTER approach, our pipeline is well 
parallelisable in the sense that each MC cycle can be run in its own 
CPU. This allows for a linear scaling with the number of CPUs 
and with the number of MC cycles performed. In our pipeline the 
most CPU time consuming step is the SDE noise generation. In an 
AIX IBM pSeries690 machine with Power4 processor running at 
clock speed of 1.1 GHz it takes ~ 1 h on a single processor to pro- 
duce one noise TOD for a single detector with length corresponding 
to 7 months of mission time (~ 2 x 10 9 samples). We have also 
implemented a parallel version of the SDE code which brings the 
required time to generate one noise TOD down to ~ 9 min with 14 
processors. 

Destriping and map coaddition of one 7 month TOD will take 
~ 50 s of CPU time in a single processor serial job. 



7 SIMULATIONS AND RESULTS 

To test our joint destriping and MASTER approach in the estima- 
tion of the angular power spectra we considered simulations rel- 
evant for the 100 GHz channels originally planned for the LFI 
instrument on-board the PLANCK satellite. The simulations in- 



8 T. Poutanen, D. Maino, H. Kurki-Suonio, E. Keihdnen, E. Hivon 



Figure 8. Top panel: Output map without destriping. Bottom panel: Out- 
put map after destriping. 



Figure 9. Top panel: Noise-only output map without destriping. Bottom 
panel: Noise-only output map after destriping. 



eluded CMB sky and instrument noise realizations. The applied 
signal+noise simulation procedure was described in Section|6| 

The scanning strategy, antenna beam response and the instru- 
ment noise characteristics applied in the simulations were those de- 
scribed in Sections l3.2l and [5] When the noise TODs were gener- 
ated the effect of using the full set of 24 radiometers was modelled 
as described in Section l7l2l When the signal TODs were generated 
a single detector scanning was assumed. 

The angular power spectrum C\ of the ACDM cosmological 
model was used as the input spectrum in the signal+noise simula- 
tions. The applied input spectrum is depicted in Fig.|3| 

Fig. IH shows a map produced from a simulated CMB+noise 
TOD, both before and after destriping. Fig.|5|shows the same for a 
pure noise TOD. 



7.1 Power Spectrum Estimate 

The power spectrum estimates Ci were solved from Eq. < 1 5 i . In a 
real CMB experiment we will not have the true angular power spec- 
trum Cf 1 available for generating the signal bias estimate S^ 10 . In- 
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Figure 10. Top panel shows the noisy power spectrum estimate (solid 
curve) used as an input spectrum to the signal-only MC to determine the 
estimate for the signal bias. The true Cf 1 is shown as well (ACDM, dashed 
curve). Bottom panel shows the obtained signal bias estimate. The sky cov- 
erage is nearly full sky. 
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Figure 11. CMB power spectrum estimated by using the destriping tech- 
nique and MASTER approach for the 100 GHz channels of PLANCK LFI. 
The use of 24 radiometers has been modelled by dividing the rms white 
noise level of a single detector with a factor \/24. Nearly full sky has been 
observed. Top panel shows the power spectra obtained for a single sky and 
noise realization: pseudo power spectrum Ct - grey points, the estimated 
power spectrum of the noise - dashed line, the estimate of the CMB power 
spectrum Ci - black points, and the input power spectrum of the true sky 
Cj h - solid line. Bottom panel shows the estimate (midpoints of the error 
bars) of the CMB power spectrum in bins of 10. The error bars are derived 
from = 450 MC realizations of signal and noise. The solid line is 

the input power spectrum. The magnitudes of the error bars are ±lo\ For 
clarity, only the odd bins are shown. 



stead a noisy power spectrum estimate can be utilized as described 
in Section |S| The power spectrum estimate and the signal bias for 
nearly full sky are shown in Fig. llQI The signal bias for the galactic 
cut can be obtained similarly. The applied noise bias (Ne) is de- 
picted in the bottom panel of Fig. Q (nearly full sky). The kernel 
matrices (M u i) derived in Section|5|were applied. 

The resulting angular power spectrum estimates and the cor- 
responding ±1<t error bars are shown in Figs. llll and ll2l for the 
nearly full sky and for the cut sky cases, respectively. Individual Ci 
are dominated by the instrument noise for I > 1000. 

The bottom panels of Figs. [TT1 and [T2l show the estimates 



E 



2nA£ 1 



(16) 



of the CMB angular power spectrum in bins b of A£ = 10. The 
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Figure 12. Same as Fig. II II except that the galactic region with |b| ^ 20° 
has been cut out from the observations. 



error bars (bottom panels) were determined from the MC variance 
of the binned power spectrum estimates. The good match between 
the estimated and the input power spectra shows that the combined 
destriping and MASTER approach produces power spectrum esti- 
mates with good accuracy. 

In order to examine the possible bias in our angular power 
spectrum estimates the mean (Ct)MC of the Cb spectra available 
from our MC realizations was determined. The relative differences 
((C 6 )mc - are shown in Figs.QJandQll 

The signal bias and the power spectrum estimates will be cor- 
related if the same seed values are used when the random sky re- 
alizations are generated in the signal-only and in the signal+noise 
simulations. This could lead to too optimistic results. We inten- 
tionally used different seed values between the signal-only and the 
signal+noise simulations. 

It is evident from Figs. ll3l and ll4l that no systematic bias can 
be observed and the variation from bin-to-bin is mainly caused by 
the residual cosmic variance and the instrument noise that remain 
after the averaging. If any bias exists its relative level is at most of 
the order of 0. 1 per cent for i < 800. 

If we had not accounted for the signal bias our method would 
have produced biased estimates. To see how large is the effect we 
did the same analysis, keeping S( = 0. This case is shown in 
Figs. \13\ and [l4| as well. A clear bias at high i (t > 1000) can 
be detected. 
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Figure 13. In order to reveal the possible bias in our method, we show 
the relative differences (thicker step curve) between the mean binned power 

450 MC realizations of signal 



spectrum estimate (Cj)mc (from A^J 11 ' 
and noise) and the binned input power spectrum Cf representing the under- 
lying theoretical spectrum. Theyaxis quantity is ((Cb)MC~ )/Cj . The 
results of the nearly full sky are shown. The symmetric smooth curves show 
the ±cr limits, where cr = A rof C b /(C*' 1 (A r ^J n) ) 1 / 2 ) is the "expected" 
standard deviation from MC noise (the reference value from Eq. I19h . No 
obvious bias is seen. The thinner step curve shows the situation without ac- 
counting for the signal bias. In this case there is a clear positive bias for 
larger £. 



7.2 Covariance Matrix and Statistics 

The covariance matrix of the power spectrum estimates was de- 
termined from the Amc = 450 signal+noise MC realizations. 
The elements of the matrix are C«/ = {(C'e — (C^)mc)(CV ~ 
(CV)mc))mc- The diagonals Ci,i+i . . . Ct,e+6 from the upper 
right triangle are shown in Fig. ll5l (galactic cut). Before plotting, 
the covariance matrix was normalized, 



.nnorm 



(17) 



to make all the main diagonal elements equal to one. The ^-to-^ 
variation of the off-diagonals is mainly caused by the MC noise. 

The second and the fourth diagonals are mostly negative. 
Their mean (up to £ = 800) values are -0.116 and -0.0265, re- 
spectively. No level shift can be detected in the sixth or higher even 
diagonals. Excluding the diagonals Cej±2 and Ce,e±4 the elements 
of the rest of the covariance matrix were buried under the MC noise. 
An analytical model for the covarianc e matrix C«/ of the pseudo- 
Ce estimators exists lEf stathio J2003I) . In the CMB dominated re- 
gion (below i ~ 800) it forecasts for the mean levels of the sec- 
ond and the fourth diagonals -0.1 13 and -0.0266, respectively. The 
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Figure 14. Same as Fig.fHlbut now the galactic region is cut out from the 
output maps. 

good match between the model values and the values obtained from 
the MC realizations is notable. The same model gives mean values 
around —1 x 1CF 4 for the rest of the diagonals shown in Fig. 1151 
No level shifts that would stand out from the MC noise could be 
detected in the covariance matrix of the nearly full sky. 

The theoretical variance for the full-sky case is given by 
Eq. For the cut sky less information is available, which in- 
creases the variance. The cut sky also introduces correlations be- 
tween the estimates Ci. The binned values Cb should be less cor- 
related. According to a "rule of thumb" (for bins Cb) the vari- 
ance should be larger approx i mately by a factor l// 3 k v (e.g. 
IScott. Srednicki & Whitelll994 iHobson & Magueiid [l996). Thus 
we shall use as reference values for the variances 



A rc fCV 



and 



cT + 



(2£ + l)f sky 



A 2 if 1 ^ f £(£ + l)A ieC C e \ 

A - fCb =zFM — a? — ) 

eeb v 7 



(18) 



(19) 



A 2 MC d e /A 2 ci d e and A 2 MC C b /A 2 c! C b , 



We show the ratios of the MC variances from our simulations 
to these reference values, 
in Fig.fT6l 

The top panel of Fig. ll6l shows that the MC variance of the un- 
binned power spectrum estimate is larger than its reference value. 
Lowering the sky coverage will produce a larger ratio. The ra- 
tio remains relatively flat up to £ ~ 1200. At higher I a posi- 
tive bump appears. The horizontal solid lines indicate the expected 



Figure 15. The diagonals jj+i . . . Ci e+6 of the covariance matrix are 
shown from top down. The matrix is normalized to C/j = 1. The expected 
std of the MC noise of the off-diagonals is (Nmc)' 1 ^ 2 = 0.047, where 
Nmc = 450. The std of the diagonals shown in the figure vary in the range 
0.046 . . .0.050. Note that the levels of the diagonals Ci^ + 2 and 
are systematically negative. The values of Cf seem to increase with 
increasing t. The sky coverage is galactic cut. 



levels derived from the analytical model of the covariance matrix 
C U i <Efstathioul2003ll The model is valid in the CMB dominated 
region. In the nearly full sky coverage the MC variances follow 
closely the model variances. In the galactic cut the model is valid 
only at high enough £ (above £ ~ 5). The MC variances of the 
galactic cut are in average slightly larger than the model variances. 
The mean relative difference is 1.8 per cent. 

The ratios between the MC variances of the binned estimates 
and their reference values are shown in the bottom panel of Fig. 1161 
Since the estimates at £ and £ + 2 are anticorrelated, binning them 
into the same bin reduces the variance by more than Al. The less 
correlated the bins are, the closer to the reference value (Eq. I19i ) 
we expect their variance to be. A smaller sky coverage implies more 
correlation and thus a higher excess variance. For our cut sky the 
remaining excess variance is 10 ... 20 per cent for bins of Al = 10 
up to £ ~ 1000. For nearly full sky no excess variance is seen below 
£ = 800. For larger I there is more excess variance. This happens 
when the spectrum becomes noise dominated. 

We determined also the statistical distributions of the angular 
power spectrum estimates obtained with our Ci estimation tech- 
nique. Examples of the PDFs (Probability Density Function) of the 
unbinned (Ci) and the binned (Cb) estimates are shown in Figs. 1171 
and !18l These PDFs were compared to the PDFs of a central xl 
model (also shown). The match between the MC PDFs and the 
model PDFs is good in the case of nearly full sky coverage. In 
the galactic cut the match is poorer especially in the PDFs of the 
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Figure 16. The ratios of the variances A^qC( / A^ cf C'i (top panel) and 
A^ lc C(,/A^ cf Cb (bottom panel). A|, c refers to the MC variance and 
A 2 £ is the reference value defined in Eqs. 1181 and 1191 . In Eq. 1181 the 

quantity (JV{)mc obtained from iVj^, = 100 noise-only MC realizations 
was used instead of (JVg). The number of signal+noise MC realizations 
was 450 when determining the variances. The horizontal lines in the top 
panel were derived from the analytic al model of the covariance matrix of 
the power spectrum estimates lEfstathiou 2003). 



Figure 17. PDFs P(x) of Cg for some selected multipoles t. The bars 
are histograms derived from 450 signal+noise MC realizations. The solid 
curves are the PDFs from the central model with v = (2( + l)/ s ky 
The std of the model is set to A re fC^, which is the reference value of the 
std of Ci defined in Eq. 1181 . In the plot the x axis quantity of the model 
PDF is first normalized with A re fQ and then the PDF is shifted along the 
x axis so that its expectation value will coincide with x = 0. 



unbinned estimates. This is as expected, since the degrees of free- 
dom (v) in the model PDFs were determined from the formula 
v — (26 + l)/ s ky, which is just a "rule of thumb" approximation 
in the galactic cut case, and does not correctly take into account the 
correlation between nearby Cg. The obtained PDFs confirm the ear- 
lier finding that the deviation of the MC variance from its reference 
value becomes larger when the sky coverage fraction is lowered 
(see Fig.ll6>. 



8 CONCLUSION 

We have demonstrated that the combination of destriping technique 
and our MASTER approach can tackle the extraction of the CMB 
angular power spectra. The approach was found to work well yield- 
ing accurate estimates of the true power spectrum of the sky. As a 
practical example, we considered the 100 GHz channels of the LFI 
instrument of the PLANCK satellite. 

Destriping as a map-making method is general in a sense that 
it requires no prior information either on the instrument noise prop- 
erties or on the actual beam shape. We can expect that even in the 
case of non-stationary noise we are able to recover the baseline 
magnitudes accurately. The MASTER technique requires knowl- 
edge of the noise characteristics. In realistic CMB experiments the 
noise parameters need to be estimated from the measured data. 



However, we did not address the noise estimation techniques in 
this study. Thus the correct instrument noise model was assumed 
throughout this paper. 

We found that the effect of the sky signal on destriping does 
not cause a significant distortion on the power spectrum estimates. 
Therefore no filter function is needed to compensate for it. 

Instead we discovered an effect related to the clustering of the 
detector pointings. Because of the high stability of the motion of 
the PLANCK satellite the centre points of the measurements taken 
during the 1 h period between satellite spin-axis repointings form 
tight clusters of 60. In our idealized simulation these were all as- 
sumed to fall on a single point, which exaggerates the effect some- 
what. This leads to a noticeable effect at high t (for I > 700 in our 
simulations). We were able to account for this effect, and remove 
it, by introducing the concept of signal bias in the Cg estimate, and 
determining it by signal-only MC simulations. We also derived a 
reasonably accurate analytical estimate for the signal bias. 

It was shown that the signal bias will be reduced when data 
from several detectors are combined, since that increases the den- 
sity of detector pointings. In additio n, due to satel lite spin axis nu- 
tation (amplitudes up to 1.5 arcmin. lvan Leeuwen et all2002l) and 
spin rate variation the detector pointings from successive scan cir- 
cles will not fall exactly on top of each other, in reality. Therefore, 
in the case of the real PLANCK experiment this signal bias may be 
small enough to be ignored. 

An ensemble of angular power spectrum estimates were pro- 



12 T. Poutanen, D. Maino, H. Kurki-Suonio, E. Keihanen, E. Hivon 



NEARLY FULL SKY 



GALACTIC CUT 



0.4 - 

0,3 
?0.2 
0.1 




-2 2 



-4 -2 2 4 



0.4 r 

0.3 
J 0.2 
0.1 



-4 -2 2 4 





-4 -2 2 4 



Figure 18. Same as Fig.rnlbut now the PDFs are for the binned estimates 
C;,. In the central xt model the number of degrees of freedom is set to 
v = ^2i£b tyfsky The std of the model is set to A rc fC;,. The bin 

size is A£ = 10. 



duced by applying our estimation method to the maps of the sig- 
nal+noise MC realizations. Our results showed that, after remov- 
ing the signal and noise biases, any remaining bias was buried un- 
der the MC variations; if any bias exists its relative level was not 
higher than 0. 1 per cent. 

There was a good match between the ±lcr error bars obtained 
from the MC realizations an d derived from t he analytical model 
of the pseudo-CV estimators lEfstathiou 2003). This indicates that 
the implementation related performance losses of our method are 
small. 
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APPENDIX A: SIGNAL BIAS 

In our signal-only simulations (Section |5J we noted an effect (at 
high I) in the angular power spectra of the output maps, that was 
not explained by beam size, pixel size, or window function, and 
could not properly be modelled in terms of a filter function. We 
studied this effect in detail, and conclude that 

(i) The effect can be modelled as an offset St, which we call 
signal bias. 

(ii) The effect is due to the resolution by which the sky is sam- 
pled by the detector pointings. 



CMB Power Spectrum Estimation with Destriping 1 3 



Since in the realistic PLANCK experiment the effect will be less 
important than in our idealized simulation, we have relegated the 
detailed discussion of this effect in this Appendix. 

We pick up the discussion from Section [5] where we intro- 
duced the signal bias in Eq. J14> . reproduced here: 



(Ct) 



S e + {N e ) . 



(Al) 



The mean pseudo spectra ((CV)mc) obtained in Section [5] were 
inserted in Eq. <A1> and the deconvolved signal bias estimate 
J2 e > M~}Sf°/Bj was solved (with {N e } = 0). It is shown in 
Fig-Elfor the ACDM model and for the nearly full sky. Fig.lA2l 
depicts similar things for the OCDM case. The galactic cut cases 
look similar but somewhat noisier. The signal bias estimates have 
a remarkable resemblance even though the cosmological models 
are widely different. This indicates that the signal bias approach is 
much better suited to our map-making method than the filter func- 
tion. 



Al Pointing Distribution 

Let us consider CMB only (no noise) and assume no baseline re- 
moval for the time being. Let k index the pixels of the final output 
map. As a result of averaging the observations the output map pixel 
temperature Tk of pixel k is 



1 

= — > Tn 
Pk ^ 



(A2) 



where pk is the number of hits in pixel k, i indexes a hit in pixel k 
and Tik is the observed sky temperature of hit i. The pixel temper- 
ature Tk can be written as 



T k = T k + T k -T k = T k + AT k 



(A3) 



where T k is the beam and pixel smoothed sky temperature of pixel 
k. The temperatures Tk and Tk can be given in terms of the expan- 
sion coefficients ai m (see Eq. 0) 



at m BiYe m (iik) 



and 

Tk z 



(A4) 



(A5) 



where B' e contains the beam and pixel smoothing of the sky map 
and Be contains the same for the output map. rife is a unit vector 
pointing to the center of the pixel k. rijfe is a unit vector of the point- 
ing of the observation i in pixel k. Performing the reverse spherical 
transform to both sides of Eq. <A3> one obtains a relation for the 
expansion coefficients 



f2 n 



fe 



T k Y e * m (rife) = Q p TkY e * m (nk) + 



fe 



+n p ) ATfel7 m (rife), 



fe 



(A6) 



where only the hit pixels will be included in the sums. 

The ensemble mean of the pseudo power spectrum Ci can be 
obtained as 



(Ci) 



•2f 



i-y Y, <K^W m (nfe) 



(A7) 
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Figure Al. Top panel: Mean deconvolved pseudo power spectrum 
M~}{Cei)wic/Bj (dashed curve), ACDM input spectrum C\ h 
(solid upper curve) and deconvolved signal bias estimate (difference of 
those two, solid lower curve). Destriping has been applied in the map- 
making stage. Middle panel: Deconvolved signal bias estimate shown with 
expanded y-scale. Bottom panel: Deconvolved signal bias estimate when 
no baseline removal has been performed (no destriping). Also shown are 
the deconvolved analytical approximation, Eq. IA22I . of the signal bias 
Sp PP (black curve) and the reference value A re f Ce divided by (N^L) 1 / 2 



(dashed curve), where Ny^-, 
produce (Ci)mc- 



MC' 

450 is the number of MC realizations to 



As a consequence of Eq. <A6> the mean can be split in three terms 

<C,) = (C?) + (Cf > + <C? A ) (A8) 
with the following formulas 



(CI) 



n 



2i + 



(C?) = £ ^<ATfeAr;)F/ m (nfe)F, m (n p )(A10) 



fl J2 ^<T fc r;)y; m (nfe)^ m (n p ) , (A9) 
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Figure A2. Same as Fig. lAlK two upper panels) but now the cosmological 
model is OCDM. 



and 



<C?A) = 27TT E E «r fc Ar;)F; m (n fc )^ m (n P ) 



m— — Z k,p 



+ (T:AT p )Y em (n k )Y e * m (n p ))(AU) 

The double sum over pixel indices k and p extends only across the 
hit pixels of the output map. Inserting Eq. <A4> and applying Eq. 
one obtains for {Cf } 

(Cf)=J2 M w B% Cf , (A 1 2) 

t' 

where M tl j is the mode-mode coupling kernel 

Mu ' = 2lTT E |fi P ^F/ m (n fe )W(n fc )| 2 .(A13) 

m,m' — — — t 1 k 

The pseudo power spectrum {Ci) can thus be expressed as a sum 
of Eq. <A12> and two additional terms 



(C e ) = M u ,BlGf + (Cf ) + {Cf A ), 



(A14) 



where we can identify Se = (Cf } + {Cf A ) as the signal bias 
defined in Eq. <14> . 



A2 Analytical Approximation to the Signal Bias 

A2. 1 Single Detector 

With the aid of some simplifying assumptions a reasonably accu- 
rate analytical estimate for St. = (Cf } + {Cf A ) (see Eqs. <A10t 



and IA1 can be derived. We consider here the case of a single 
detector and nearly full sky coverage. 

Let us consider an approximation to (C A ) first. This is called 
S A . There exists a non-zero correlation between the temperatures 
Tk of different output map pixels. The temperature ATk (see Eqs. 
<A2t and <A3» is a linear combination of differences between the 
temperatures of the hits falling into pixel k and the temperature at 
the centre point of that pixel. The magnitudes of the pixel-to-pixel 
correlations are smaller for these differences than for the tempera- 
tures themselves. Thus we can apply the following approximation 



(ATfcAT;) = {\AT k \ 2 )5 kp . 



(A15) 



By applying the summation rules of the spherical harmonics we 
obtain for Sf 



ST 



(A16) 



The summation extends over the hit pixels of the output map. The 
number of those pixels is n P i X . The formula exhibits no £ depen- 
dence. The expectation value (jATfe| 2 ) can be evaluated with the 
help of the input spectrum Cj h 



(A17) 



where 



x k _ B e B e 



P,(cos(4)) 



IBgB'i 
p k 



^ft(cos(^)) + J B, 



(A 18) 



The upper limit for I is £ m ax, Pl{x) is the Legendre polynomial, 
cos(0y) = n lk ■ n jk and cos(6» l fc ) = n ik ■ n k . 

By applying the pointings (n^) of the 100 GHz LFI detector 
scanning (see Section[5} S^ was calculated and the result is shown 
in Fig.[A3]for ACDM and OCDM. 

Let us examine the approximation of {Cf A ) next. It is denoted 
by Sf A . We shall insert Eqs. JA2), JA4} and jA5j in Eq. iATTl 
and use an approximation that is applicable in the nearly full sky 
coverage 



v / J Ytmjii^Ytt m i (rtfc) ~ 5 ee i5„ 



(A19) 



(A20) 



The summation is extended across the hit pixels. We obtain for Sf 

Sf A — 2C) h (BiB' e Ue — -B?) , 
where 

Pk 

(A21) 



t^pix P k 

k i=l 



Finally we can write the analytical approximation (S^ pp ) of 
the signal bias Si, 

S* pp = Sf + Sf A . (A22) 

It is shown in Fig.[A3]for both ACDM and OCDM. For low i (£ < 
400) the signal bias is dominated by Sf and for high £ by S A . 
At low multipoles the level of the signal bias is several orders of 
magnitude lower than the level of the CMB power spectrum. At 
high £ their levels become closer to each other. 
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Figure A3. Analytical approximations of the total signal bias S^ pp due to 
the detector pointing distribution (solid lower curves) and the ^-independent 
part of it (dashed curves). Top panel is for ACDM and bottom panel 
for OCDM. The respective beam and pixel weighted CMB angular power 
spectra are also shown (solid upper curves). 



As shown in the bottom panel of Fig. lAll there is a good match 
between the analytical approximation of the signal bias S^ pp and 
the signal bias estimate Sf 10 derived from the MC realizations. 

One might have expected that increasing the map resolution 
(decreasing the pixel size) would have reduced the high-£ tail of the 
signal bias (and the filter function) or at least would have shifted it 
to higher multipoles. Contrary to this expectation the high-^ tail re- 
mains almost unaltered when the map resolution is increased (see 
Fig. |A4j. We can now explain this. The temperature differences 
Tik — Tfe in average and the pixel size Sl p become smaller when 
the pixel size is reduced. This would imply decreased ATt and 
S^. However, when the pixel size is reduced there will be in aver- 
age lower number of hits per pixel and those hits will be more un- 
evenly spread inside the pixel area. An uneven spread would yield 
increased ATI-. The net effect of these two opposite phenomena is 
that the magnitude of remains nearly unaltered in this case and 
the high-£ tail remains roughly in the same place. 

In other words, the crucial resolution element here is not the 
pixel size, but the distance between neighboring detector pointings, 
separated by 3.3 arcmin along the ring, the rings themselves being 
separated by 2.5 arcmin. 



A2.2 Multiple Detectors 

Destriping can be extended straightforwardly from a single detec- 
tor TOD to multiple TODs produced by several detectors. The base- 
lines can be determined jointly for all detectors and subtracted from 
the TODs. A single multidetector map can be obtained by consider- 
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Figure A4. Same as the bottom panel of Fig. IA II but now the input map 
and output map had higher resolutions (iV s j^ e = 2048 and N s \^ c = 1024, 
respectively). 

ing the hits of all detectors in a pixel. The pseudo power spectrum 
can be extracted from the obtained multidetector map. 

The signal bias will be altered when data from several detec- 
tors are combined, since that increases the density of detector point- 
ings. To study this issue the pointing matrices of 70 GHz detectors 
of PLANCK LFI instrument were available. There are totally six 
LFI detectors at 70 GHz. Each single detector pointing matrix cov- 
ered 12 months of observations. We calculated the analytical ap- 
proximation Sf pp of the signal bias by utilizing Eqs. <A16t . IA17I . 
<A18L <A20t and <A2H . The number of detectors was modelled by 
extending the calculations over the respective number of pointing 
matrices. It should be noted that due to differences in the sky cover- 
age, in the mission time and in the sampling rates the absolute levels 
of the signal biases of LFI 70 GHz and 100 GHz detectors cannot 
be directly compared here. The results obtained with 70 GHz detec- 
tors reveal the relative impact that the number of detectors will have 
on the signal bias of the pseudo power spectrum extracted from a 
multidetector map. 

The signal biases of one, three and six detectors are depicted 
in Fig. lA5l At low £ (£ < 300) the signal bias is nearly independent 
of the number of detectors. The absolute level of the signal bias 
is highest at this range of multipoles. At high £ (£ > 1200) the 
signal bias is clearly dependent on the number of detectors. At this 
range of multipoles the signal bias is dominated by Sjr' . The ratio 
between one detector and three detector Sf" is 3.15. The respective 
ratio between one detector and six detector is 5.83. Thus the 
high-^ level of Sg PP is nearly inversely proportional to the number 
of detectors involved. At the medium range of I's the behavior is 
more complex due to the abrupt change of sign of Sf A . 

In light of this information, the level of the signal bias of a 
LFI 100 GHz detector can be compared to the typical level of the 
angular power spectrum of the instrument noise. The absolute level 
of the signal bias is highest at low £ and this part is nearly indepen- 
dent of the number of detectors involved. Comparing the analytical 
signal bias of Fig. I A3 1 (single detector case) to the angular spec- 
trum of the instrument noise (shown for the modelled full set of 24 
detectors in the bottom panel of Fig.0 one can notice that in this 
case the low-^ level of the noise bias is at least couple of orders of 
magnitude higher than the respective level of the signal bias. 

A3 Estimating the Signal Bias 

In a real CMB experiment we will not have the true angular power 
spectrum C\ h available for estimating Si. Since we have concluded 
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Figure A5. Analytical approximations (S^ pp ) of the signal biases of one 
(A), three (B) and six (C) combined detector signals. The bias of one detec- 
tor is the mean of the signal biases of the individual 70 GHz LFI detectors. 
The bias of three detectors is obtained by considering the pointings of 70 
GHz LFI detectors 18, 19 and 20. The pixel size of the output map was ac- 
cording to A r s j ( j c = 512 and all the detector pointings fell in the centres of 
N s idc = 1024 pixels. The horizontal lines show the contribution of in 
each case. Dashed lines indicate negative values of S^ pp . The cosmological 
model was ACDM. 

that its dependence on C| is weak, this is not a concern. In place 
of G\ h we can use the noisy power spectrum estimate Ce that is 
obtained from the real data. 

In order to produce the noisy power spectrum estimate a single 
realization of the signal+noise pseudo power spectrum Ci was gen- 
erated by following the procedure described in Section|6| A single 
100 GHz LFI detector was considered. The estimate Ci was solved 
from Eq. J 1 5 1 by setting St = and by using (Ni) determined ac- 
cording to the pure noise MC procedure (see Section|6j. The ob- 
tained C e was binned with A£ = 10 bins at I ^ 1200 and with A£ 
= 50 bins at I > 1200 and smoothed by a cubic spline (Fie. 1101 top 
panel). 

The noisy power spectrum estimate was used as the input 
spectrum for signal-only MC simulations (N^^ — 450) to pro- 
duce the estimate of the signal bias (Fig. 1101 bottom panel). De- 
striping was applied in the map-making. As expected, the signal 
bias has a remarkable resemblance to the signal bias obtained from 
simulations with the true ACDM spectrum (see Fig. lAH . 
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